The localization transition at finite temperatures: electric and thermal transport 
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The Anderson localization transition is considered at finite temperatures. This includes the elec- 
trical conductivity as well as the electronic thermal conductivity and the thermoelectric coefficients. 
An interesting critical behavior of the latter is found. A method for characterizing the conductivity 
critical exponent, an important signature of the transition, using the conductivity and thermopower 
measurements, is outlined. 
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I. INTRODUCTION 

Anderson localizationi'2. is a remarkable, and very 
early, example of a quantum phase transition (QPT), 
where the nature of a system at T = changes abruptly 
and nonanalytically at a point, as a function of a con- 
trol parameter. Here it implies that the relevant quan- 
tum states of the system acquire an exponential decay 
at large distances (similar to, but much more complex 
than, the formation of a bound state). In the original 
paper this phenomenon was discovered via a change in 
the convergence properties of the "locator" expansion. 
In the disordered tight-binding model, with short-range 
hopping, this is the expansion about the bound atomic, 
or Wannier-type orbitals. As long as the expansion con- 
verges, the relevant eigenstates are localized; its diver- 
gence signifies the transition to "extended", delocalized 
states. This was reviewed, hopefully pedagogically, in 
Ref. 0]. Later, Mott introduced the very useful picture 
of the "mobility edge" within the band of allowed en- 
ergies, separating localized from extended states^^i In the 
lower part of the band, the states below the "lower mobil- 
ity edge" are localized, while those above it are extended. 
When the disorder and/or the position of the Fermi level, 
Ep, are changed, a point where they cross each other is 
where the states at Ep change nature from extended to 
localized, and this is a simple and instructive model for 
a metal-insulator transition at vanishing temperature T. 

The analysis of the above localization transition cen- 
ters on the behavior of iTf){E), the conductivity at energy 
E, which would be the T = conductivity of the sample 
with Ep — E. ao{E) vanishes for E below the lower mo- 
bility edge Em and goes to zero when E approaches Em 
from above^ 



aoiE)=A{E-E„ 



(1) 



The characteristic exponent for this, x, is an important 
parameter of the theory. It is expected to be universal 
for a large class of noninteracting models, but its value 
is not really known, in spite of the several analytical, 
numerical and experimental methods used to attempt its 
evaluation. 

The electron-electron interaction is certainly relevant 



near this transition, which may bring it to a different uni- 
versality class. This is a difficult problem. The bench- 
mark treatment is the one by Finkel'stein.^ Since it is 
likely that the situation in at least most experimental 
systems is within this class, it would appear impossible 
to determine the value of x for the pure Anderson tran- 
sition (without interactions))^ 

In this paper we still analyze the thermal and ther- 
moelectric transport for a general model with ctq^E) be- 
having as in Eq. ([T|). This is certainly valid for nonin- 
teracting electrons and should be valid also including the 
interactions, as long as some kind of Landau Fermi-liquid 
quasiparticles exist. In that case (Tq{E) for quasiparticles 
is definable and thermal averaging with Fermi statistics 
holds. Even then, however, various parameter renormal- 
izations and interaction corrections^ should come in. The 
effect of the interactions on the thermopower for a small 
system in the Coulomb blockade picture was considered 
in Ref. Q, and correlations were included in Ref. 
Even in the latter case, the simple Cutler-MottiS' for- 
mula (Eq. [221 derived in section UlI)) was found to work 
surprisingly well. 

It should be mentioned that the sharp and asymmetric 
energy-dependence of ao{E) near the mobility edg o^°i^^ 
should and doe a^°i^^'^'^ lead to rather large values of the 
thermopower. Exceptions will be mentioned and briefly 
discussed later. Large thermopowers are important for 
energy conversion and refrigeration applicationsi^ and 
this clearly deserves further studies. 

A serious limitation on the considerations presented 
here is that the temperature should be low enough so that 
all the inelastic scattering (electron-phonon, electron- 
electron, etc.) is negligible. For simplicity we consider 
here only longitudinal transport (currents parallel to the 
driving fields). Thus, no Hall or Nernst-Ettinghausen 
effect! We also do not consider the thermoelectric trans- 
port in the hopping regime. It might be relevant even in 
the metallic regime (chemical potential fi > Em), once 
ksT > fi — Em- This may well place limitations on the 
high temperature analysis we make in the following. 

In section |TT] we review the basic concepts behind the 
scaling theory^ for the transition, and reiterate the criti- 
cal behavior of the conductivity as in Eq. ([T]), obtaining 
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also its temperature dependence. In section Hill we derive 
all results for the thermal and thermoelectric transport 
and analyze the scaling critical behavior of the latter as 
function of temperature and distance from the transition. 
A brief comparison with experiment is done in section 
IIVI and concluding remarks are given in section [V] In 
the appendix we present a proof that the heat carried 
by a quasiparticle is equal to its energy measured from 
the chemical potential, /i, hoping that this also clarifies 
the physics of this result. This derivation is valid for 
Bose quasiparticles (phonons, photons, magnons, exci- 
tons, plasmons, etc.) as well. 

II. THE ZERO AND FINITE TEMPERATURE 
MACROSCOPIC CONDUCTIVITY AROUND 
THE ANDERSON LOCALIZATION TRANSITION 

A. The Thouless picture within the 
tunnel-junction model 

We start this section by briefly reviewing the tunnel- 
junction picture of conductio n^^i^^ at T = 0, which is a 
useful way to understand the important Thoulessii pic- 
ture for such transport. Consider first two pieces (later 
referred to as "blocks") of a conducting material with a 
linear size L, connected through a layer of insulator (usu- 
ally an oxide) which is thin enough to allow for electron 
tunneling. The interfaces are assumed rough, so there is 
no conservation of the transverse momentum: each state 
on the left interacts with each state on the right with a 
matrix element t with a roughly uniform absolute value. 
The lifetime for an electron on one such block for a 
transition to the other one is given by the Fermi golden 
rule (at least when tunneling is a weak perturbation): 

r-' = ^WNr{E,), (2) 

where |ip is the average of the tunneling matrix element 
squared and Nr{Ep) is the density of states on the final 
(right-hand) side. Taking the density of states (DOS) 
in the initial side to be Ni{Ep), we find that when a 
voltage V is applied, eVNi{Ep) states are available, each 
decaying to the right with a time constant r^^ , so that the 
current is 

I = e^m{E^)T;'V, (3) 
and the conductance is 

G = e'NeiE,)/T, = ^t^N<,iE,)NriE,) , (4) 

which is an extremely useful result. This equality is well- 
known in the tunnel junction theory J^ii^ Clearly, G is 
symmetric upon exchanging I and r, as it should.— Note 
that Eqs. ^ and (|H) are valid in any number of dimen- 
sions. An important remark is that Eq. ^ necessitates 
a continuum of final states, while the final (RHS) block 



is finite and has a discrete spectrum. One must make 
the assumption that the interaction of that system with 
the outside world leads to a level broadening larger than, 
or on the same order of, the level spacingii^r^ This is 
the case in most mesoscopic systems. One then naively 
assumes that this condition converts the spectrum to an 
effectively continuous one (the situation may actually be 
more subtle). Otherwise, when levels really become dis- 
crete, one gets into the really microscopic (molecular) 
level. 

The result of Eq. ^ is very general. Let us use it 
for the following scaling picturSiii: Divide a large sample 
to (hyper) cubes or "blocks" of side L. We consider the 
case L ^ £,a; £ being the elastic mean free path and a 
the microscopic length. The typical level separation for a 
block at the relevant energy (say, the Fermi level), d^^, is 
given by the inverse of the density of states (per unit en- 
ergy) for size L, Nl{Ep). Defining an energy associated 
with the transfer of electrons between two such adjacent 
systems by Vl = ttS/t^ (r^ is the lifetime of an elec- 
tron on one side against transition to the other side) the 
dimensionless interblock conductance = G^/{e'^/iTh) 
is: 

9.=Vjd, (5) 

i.e. is the (dimensionless) ratio of the only two rele- 
vant energies in the problem. The way Thouless argued 
for this relation is by noting that the electron's diffusion 
on the scale L is a random walk with a step L and char- 
acteristic time , thus 

Note that as long as the classical diffusion picture holds, 
Dl is independent of L and tl = LP'/D, which is the 
diffusion time across the block. It will turn out that the 
localization or quantum effects, when applicable, cause 
Dl to decrease with L. For metals, the conductivity, tr^, 
on the scale of the block size L, is given by the Einstein 
relation a = DLe^dn/ (where ^ is the chemical poten- 
tial and dn/dn is the density of states per unit volume), 
and the conductance in d dimensions is given by: 

G,^a,U'-\ (6) 

Putting these relations together and remembering that 
N^{Ef) ^ L'^dn/dn, yields Eq. dS]). 

To get some physical feeling for the energy ft./r^ we 
note again that, at least for the weak coupling case, the 
Fermi golden rule yields Eq. ([2]) or: 

K=2^WK- (7) 

Thus, is defined in terms of the interblock matrix el- 
ements. Clearly, when the blocks are of the same size, 
Eq. d?]) is also related to the order of magnitude of the 
second order perturbation theory shift of the levels in 
one block by the interaction with the other. For a given 
block this is similar to a surface effect - the shift in the 
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block levels due to changes in the boundary conditions 
on the surface of the block. Indeed, Thouless has given 
appealing physical arguments for the equivalence of Vl 
with the sensitivity of the block levels to boundary con- 
ditions. This should be valid for L much larger than £ 
and all other microscopic lengths. 

Since in this scaling picture the separations among the 
blocks are fictitious for a homogeneous system, it is clear 
that the interblock conductance is just the conductance 
of a piece whose size is of the order of L, i.e., this is the 
same order of magnitude as the conductance of the block 
itself. 

The latter can also be calculated using the Kubo lin- 
ear response expression. It has to be emphasized that the 
Kubo formulation also applies strictly only for an infinite 
system whose spectrum is continuous. For a finite sys- 
tem, it is argued again that a very small coupling of the 
electronic system to some large bath (e.g. the phonons, 
or to a large piece of conducting material) is needed to 
broaden the discrete levels into an effective continuum. 
Edwards and Thouless," using the Kubo-Greenwood for- 
mulation, made the previously discussed relationship of 
Vl with the sensitivity to boundary conditions very pre- 
cise. 

The above picture is at the basis of the finite-size 
scaling^i^i theory of localization. It can also can be 
used for numerical calculations of g{L), which is a most 
relevant physical parameter of the problem, for non- 
interacting electrons, as we shall see. Alternatively, Eq. 
([7]) as well as generalizations thereof can and have been 
used for numerical computations. Powerful numerical 
methods exist to this end.^^ 

It is important to emphasize that 3> 1 means that 
states in neighboring blocks are tightly coupled while 
(7i ^ 1 means that the states are essentially single-block 
ones. gL is therefore a good general dimensionless mea- 
sure of the strength of the coupling between two quan- 
tum systems. Thus, if 5^ ^1 for small L and 5^, — for 
L — 00, then the range of scales L where ^ I gives 
the order of magnitude of the localization length, ^. 

Although the above analysis was done specifically for 
non-interacting electrons, it is of greater generality, g^ 
(with obvious factors) may play the role of a conductance 
also when a more general entity (e.g. an electron pair) is 
transferred between the two blocks. The real limitations 
for the validity of this picture seem to be the validity of 
the Fermi-liquid picture and that no inelastic effects (e.g 
with phonons or electron-hole pairs) occur. 

The analysis by Thouless^ of the consequences of Eq. 
([5]) for a long thin wire has led to extremely important 
results. First, it showed that ID localization should man- 
ifest itself not only in "mathematically ID" systems but 
also in the conduction in realistic, finite cross-section, 
thin wires, demonstrating also the usefulness of the block- 
scaling point of view. Second, the understanding of the 
effects of finite temperatures (as well as other experimen- 
tal parameters) on the relevant scale of the conduction, 
clarifies the relationships between g(L) and experiment 



in any dimension. Third, defining and understanding 
the conductance g{L) introduces the basis for the scal- 
ing theory of the Anderson localization transition.— Here, 
we use the results for the (macroscopic) T = conduc- 
tivity around the localization transition to get the finite 
temperature conductivity there. 

B. The critical behavior of the T — conductivity 

Near, say, the lower mobility edge, Em, the conductiv- 
ity a'Q{E) vanishes for E < Em and approaches zero for 
E Em from above, in the manner: 

ao{E) ^ AiE - Emf , (8) 

A being a constant and x the conductivity critical ex- 
ponent for localization, which has so far eluded a pre- 
cise determination either theoretically or experimentally. 
Within the scaling theory,— x is equal to the critical ex- 
ponent of the characteristic length (^) , because 



In that case, an appealing intuitive argument by Mottj^i 
and Harris^^- places a lower bound on x: 

X > 2/d. (10) 

In fact, Eq. ^ may be expected to hold on dimen- 
sional grounds for any theory which does not generate 
another critical quantity with the dimension of length. 
This should be the case for models which effectively do 
not have electron-electron interactions. With electron- 
electron interactions, for example, we believe that the 
critical exponent for the characteristic length should sat- 
isfy an inequality such as Eq. (ITUl) i^ However, this may 
no longer be true for the conductivity exponent. 

C. The conductivity at finite temperatures 

In Eqs. ^ and we calculated, at T = 0, the current 
in an infinitesimal (linear response) energy strip of width 
eV around the Fermi energy. Generalizing this to an 
arbitrary energy at finite temperature, we find that the 
current due to a strip dE at energy E is 

I{E)dE = eNe{E,)T;\E)[fi{E) - fr{E)]dE, (11) 

fi{E) {fr{E)) being the Fermi function at energy E at the 
left (right). The total current is obtained by integrating 
Eq. ((TT|) over energy. For linear response fi{E) — fr{E) = 
eV[-§^]. This gives 

<yiT)^ fjEa,{E)[-^l (12) 

where cro{E) = {e'^/nh) L^'^-d.) ^-j^g conductivity 

(using Eqs. [5] and IH]) at energy E, which would be the 
T — conductivity of the sample with Ep = E. 
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D. Analysis of a{T, fi — Em) 

From now on we assume Eq. ([8]) to hold. Measuring 
all energies from the chemical potential /i and scaling 
them with T, we rewrite Eq. (I12p in the manner (we 
shall employ units in which the Boltzmann constant, ks 
is unity, and insert it in the final results) 

aiT,^i-E^)^AT-J:{[^i~Em,]/T), (13) 

where the function T,{z) is given by: 

S(z) = J dy{y + z) [ — ]. (14) 

Fig. [1] shows a numerical evaluation of this integral. 

Let us consider the low and high temperature limits 
of this expression, that will also be useful later for the 
analysis of the thermopower. 

At low temperatures, we can use the Sommefeld ex- 
pansion, to obtain: 

moUT)/A = {^l- EraT + -^T''x{x - - Emf-^. 

(15) 

Notice that is negative for a; < 1: this comes about 
since in this case the function (t{E) is concave. 

At high temperatures, one can set Em = 0, since the 
contribution to the integrals comes from energies smaller 
or of the order of the temperature, and T 3> — Em- 
Therefore we have to evaluate: 
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FIG. 1: The integral of Eq. (|14p was evaluated numerically, 
for different values of the exponent x, for — Em = 1. Shown 
is the dependence of the conductivity, measured in units of A 
(see Eq. ((SI) as a function of temperature (the energy scale is 
set by /X — Em)- At low temperatures, the conductivity satu- 
rates to a value given by Eq. (|15p . while at high temperatures 
it scales as T^, see Eq. (|16|l . The derivative at low temper- 
atures can be positive or negative, depending on x (see Eq. 



an,,n(T)IA^ E^{-^)dE. (16) 

Thus, at high temperatures (Thigh /A ^ T^, with the 
coefficient given by 2-i-cosh(q) section IIIICI we 

show that this integral can be connected with the Rie- 
mann Zeta function, and its value is given by Eq. (|36l) . 

This scaling could be used for a determination of the 
exponent x. However a much closer determination of that 
exponent would follow from the scaling of both the con- 
ductivity and the thermal and thermoelectric transport 
coefficients, which will be studied in the coming sections. 

III. THERMAL AND THERMOELECTRIC 
TRANSPORT 

A. General relationships 

Consider now the case where both a voltage V and a 
temperature difference AT are applied between the two 
blocks. We choose for convenience /c^ = 1. Both are 
small enough for linear response to hold. Here we have 
to replace the Fermi function difference in Eq. ([TT|) by 

f,E)-friE)^eV[-^]+AT[-^]. 

Then generalizing Eq. (IT^ yields for the electrical cur- 
rent 

/ = / dE^{[-e^]V + [-^]AT}, (17) 

with G given in Eq. ([6|). The first term is the ordinary 
ohmic current and the second one is the thermoelectric 
charge current due to the temperature gradient. 

Next, we derive the heat current. The heat carried by 
an electron with energy E (measured from the chemical 
potential, /i) is equal to E. This is shown, for example in 
Ref. [lO] by noting that the heat is the difference between 
the energy and the free energy. Sivan and Imrj^ii verified 
it in their Landauer-type model by calculating the flux of 
TS along the wire connecting the two reservoirs. In the 
appendix, we obtain the same result in our block model, 
from the time derivative of the entropy of each block. 
Thus, we obtain the heat current /g, 

I,^jEdE^{[^e^]V+[-^]AT}. (18) 

Here the first term is the thermoelectric heat current due 
to the voltage, while the second one is the main contri- 
bution to the usual electronic thermal conductivity n. 

In this model the ratio of thermal to electrical conduc- 
tivities is of the order of {kB/e)'^T. This is because a 
typical transport electron carries a charge e and an exci- 
tation energy of the order of fc^T and the driving forces 
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are the differences in eV and ksT. This ratio is basically 
the Wiedemann-Franz law. 

It is convenient to summarize Eqs. (fTTj) and ([T8|) in 
matrix notation^^: 



L21 L22 



V 

AT 



(19) 



where the coeiBcients Lij can be read off Eqs. (fT7|) and 



(20) 



T8|) . Since / is a function of E/T, we see that 



a/ 

dT 



Edl 

T dE 



Therefore, the two "nondiagonal" thermoelectric coeffi- 
cients: the one relating / to AT, L12, and the one relating 
Iq to V ^ L21, are equal within a factor T. 



L21/T. 



(21) 



This is an Onsager— relationship, which holds very 
generally for systems obeying time-reversal symmetry 
(and particle conservation - unitarity). The case where 
time-reversal symmetry is broken, say by a magnetic 
field, is briefly discussed in the next subsection. 

We conclude this subsection by defining and obtain- 
ing an expression for the absolute thermoelectric power 
(henceforth abbreviated as just "thermopower" ) of a ma- 
terial. Suppose we apply a temperature difference AT 
across a sample which is open circuited and therefore no 
current can flow parallel to AT. To achieve that, the 
sample will develop a (usually small) voltage V, so that 
the combined effect of both AT and V will be a van- 
ishing current. From Eqs. (IT7|) . (flQl) and (l20t we find 
that the ratio between V and AT, which is defined as 
the thermopower, S, is given by 



S 



V 

AT 



L12 



JdEEaoiE)§l 
eTjdEao{E)§l 



(22) 



B. Onsager relations in a magnetic field 

From time-reversal symmetry at H = and unitarity 
(particle conservation) follows the Onsager relation-22r— 
for the T = conductance 



(t{E,H) = aiE,-H) 



(23) 



This can be proven for our model from the basic symme- 
tries of the interblock matrix elements. This symmetry 
obviously follows for the temperature-dependent electri- 
cal and thermal conductivities a{T) and k{T). 

For the nondiagonal coefficients, the usual Onsager 
symmetry reads 



Li2{H)^L2i{-H)/T. 



(24) 



In our case, since the nondiagonal coefficients are ex- 
pressed as integrals over a symmetric function (Eq. (|23|) ). 
they also obey 



Lij{H) = Lij{-H). 



(25) 



i.e. the nondiagonal coefRcients are symmetric in H as 
well. 



C. Analysis of the thermopower 



Eq. ((22|) for the thermopower is identical to the one 
derived in two-terminal linear transport within the Lan- 
dauer formulation in Ref. [T]| . which is equal in the 
appropriate limit to the Cutler-Mott^i^Z expression: 



^,)ao{E){-m 



ea{T)T 



(26) 



where /i is the chemical potential, (To(£') is the conduc- 
tivity for carriers having energy E and a is the total 
conductivity. The Physics of this formula is clear for the 
(Onsager-dual) Peltier coefficient: a carrier at energy E 
carries an excitation energy (similar to heat, see the ap- 
pendix) of E' — /i. 

Clearly, electrons and holes contribute to S with op- 
posite signs. S will tend to vanish with electron-hole 
symmetry and will be small, as happens in many metals, 
especially in ordered ones, when the variation in energy 
of cro{E) around /z is weak. 

Having a strong energy dependence of ctq (E) , and be- 
ing very different above and below fi will cause relatively 
large values of S. We believe that this is what happens 
in disordered narrow-gap semiconductors, which feature 
in ma ny pre sent-day good thermoelectrics. As noted in 
Refs. |10llll| . the Anderson metal-insulator transition (or 
at least its vicinity) offers an almost ideal situation for 
large thermopowers. There, uq{E) vanishes below the 
mobilty edge Em (for electrons) and approaches zero, 
probably with an infinite slope, above it. Hopping pro- 
cesses in the localized phase are not considered here. A 
brief analysis in Ref. [ll| demonstrated that S scales with 
z = (/i - Em)/T: 



i» — Em . 
T 



(27) 



(Y being a universal scaling function) and assumes the 
two limits: 



and 



S ^ {H~E,„)~\ for z > 1); 



S const — z, for z ^ 1. 



(28) 



(29) 



Of course, there is no "real" divergence of Sr^ since 
when (/i — Em) (and the slope of S{T) diverges), it 
will eventually become smaller than T and the large-slope 
linear behavior will saturate as in Eq. (|29l) . 

Fig. [5] shows a numerical evaluation of Eq. (^5]) . 
demonstrating the linear low temperature regime and the 
saturation at high temperatures. Let us now make a more 
thorough investigation of the low and high temperature 
regimes. 
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For low temperatures, we can use, as before, the Som- 
merfeld expansion for the nominator and denominator, 

to obtam: ^ / (34) 



Slow 



^Tx{n - E^f-^ + 0{T^) Defining C = ^^dm, we find that: 

(30) roo 2to^-1 

A more complete expression is give in Eq. m below. C - G//3 = / _ rfm = (35) 



Thus, at T ^ /i — Em, the thermopower is linear in 
temperature: therefore: 



5, 



low 



3e(^ - £■„) 



o(r3 



(31) G(/3) = /3C(1 - 1/2^-1) = /3C(/3)r(/3)(l - l/2'5-i). (36) 



Fig. [3] compares this expression with the numerically 
evaluated slope. 

At high temperatures, one can set E,n — 0, as before. 
Therefore we have to evaluate: 



Jo 



-E^+\-§l)dE 



^ _ Jo ^ y-mi"'^ 

- eT J- EH-§l)dE 



(32) 



We can write E'^{-§^)dE = Tl^+^G{/3), with the 
dimensionless function G(/3) defined as: 



2 + cosh(m' 



(33) 



G{x+2) 



We then have Shigh - eOix+i) - 

In fact, the integral of Eq. ([55)) can be related to the 
Riemann Zeta function C: 



This gives an exact formula for the thermopower at high 
temperatures: 



high 



(l + x) 



C(l + x)(2^~l) 
eC(a;)(2^-2) 



(37) 



At a; = 0, one obtains S = 21og(2), while for x ;2> 1, one 
obtains Shigh ~ 1 + x. 

Actually, understanding the behavior for large x is 
simple: If we were to approximate the derivative of the 
Fermi-function by e~-^/^ , we would have G{(3) = r(/3), 
where F stands for the Gamma function. Then, by its 
properties, we immediately have that Shigh ~ (1 + x)/e. 

It turns out that a good approximation to Shigh (x) can 
be obtained by interpolating the exact a; = result and 
the large x result, by the form: 



s 



high 



[21og(2) 



(38) 




DO 2 



a> 

Q. 

o 
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FIG. 2: The integrals of Eq. (|26p were evaluated numerically, 
for different values of the exponent x. At low temperatures, 
the thermopower is linear with temperature, while at high 
temperatures it saturates, to a value which depends on x. 



FIG. 3: At low temperatures, the thermopower is described 
by Eq. (|30|) . The plot shows a comparison between this ex- 
pression and the evaluation of the slope extracted from the 
numerics demonstrated in Fig. [J] 
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Fig. |4] compares the exact saturation values of Eq. (p7)) 
with this approximate form. We found that the difference 
for all values of x is less than 6 percent, and therefore 
Eq. psp provides a practical working formula for the 
saturation value of the thermopower. 

An interesting feature of the crossover from the low to 
high temperature regime, is the possibility of an inflection 
point in the thermopower dependence. Similar to the 
behavior of the conductance, which grew for x > \ but 
diminished for x < 1, here there will be an inflection 
point for x <1. To see this, we have to calculate the next 
order in the Sommerfeld expansion in the nominator Q 
of Eq. ((30|). to obtain, up to corrections of order O(T^): 



Q = ^Tx{^l^E^r-^+^^T''x{x^l){x-2){^i-E^r-\ 

(39) 

This leads to the following low temperature correction 
of the thermopower: 



2 rj. 4 

Slon, « . / ^ + — x{x-l){x-7)- -— +0(r5), 



ie{n-E^) 45 



implying an inflection point for .x < 1. 



(40) 



IV. BRIEF DISCUSSION OF EXPERIMENTS 

Large thermopowers that are linear in the tempera- 
ture, at least in the metallic regime, were already found 
in the pioneering extensive work on Cerium sulfide com- 
pounds by Cutler and Leavy^i^ and analyzed by Cut- 
ler and Mott.lS It is interesting to address specifically 



the behavior around the localization transition. An ex- 
periment performed on In20^-x (both amorphous and 
crystalline), approaching the Anderson MIT, shortly af- 



confirmed qualitatively the main features 
13 Values of S exceeding 100^^ 



ter Ref. 

of Eqs. (UHl) and (gHl)-^^ Values of S exceeding ^ 
were achieved. It should be kept in mind that for a good 
determination of the critical exponent x, one needs data 
at low temperatures and small /i — E,-a . Data too far from 
the QPT, which is at both T = and /i - £',„ = 0, will 
not be in the critical region and may be sensitive to other 
effects, as will be discussed later. 

It has been customary to use only the low temperature 
conductivity to determine the critical exponent x. Using 
similar In203-x samples, the conductivity was extrapo- 
lated in Ref. (3J| to T = and those values were plotted 



against a control parameter which should be proportional 
to /i — Em when both arc small. A value of x — .75 — .8 
was found. 

It would be much better to use both the above conduc- 
tivity values and the slopes of S{T) near T — 0, according 
to Eq. (|31l) . An even better way to do that would be to 
eliminate the control parameter /x — En from Eqs. (jS]) 
and (HH]), getting 



dS 



-HT = 0)] 



-l/x 



(41) 



not having to determine the additional parameter fi — E^^ 
for each case. The data allowed us to effect this only 
approximately, see Fig. [SJ giving x = 1 ± .2. However 
getting near the QPT, this procedure appears to be the 
one of choice. 

Obviously, using the two full functions cr(/i — E„nT) 




FIG. 4: At high temperatures, the thermopower saturates, at 
a value dependent on x. The plot shows the saturation values 
extracted by evaluating Eq. (|26|) numerically (see Fig. [2]), and 
a linear dependence corresponding to Eq. (PS)) . 




log (a) 

FIG. 5: A small set of data, which should be taken as a pre- 
liminary to more extensive studies, was used to find the crit- 
ical exponent x. Eqs. ((8]) and (|28p show that the linear (low- 
temperature) regime of the thermopower has a slope which 
has an inverse power-law dependence on the distance from 
the transition. The slope of the log-log plot gives x ^ 1.1, 
with a significant error. 
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and S{fi — Ejn, T), in the critical region would place even 
more strict constraints on x. Below, we do this for the 
existing data, to demonstrate the method. Their scaling 
works well, but the value of x obtained is not likely to be 
the real critical value. This is due to a few caveats which 
will be mentioned. 

Fig. [S] compares the above predictions to the exper- 
imental data, taking the exponent x and the values of 
Em as fitting parameters. Fig. [7] shows the approximate 
data collapse obtained by rescaling the temperature axis 
of each of the measurements (corresponding to the appro- 
priate value of Em) as in Eq. (I?f)) . and the theoretical 
curve corresponding to x = 0.1. 

The fit is certainly acceptable. However, the value of 
X = .1 is both in disagreement with the previously deter- 
mined value and impossible for noninteracting electrons, 
since there x > 2/3 in three dimensions. Although, as 
explained at the end of subsection III Bl this constraint 
may not be valid with interactions, we do not take this 
last value of x seriously. Since the saturation value was 
shown be approximately given by a x-l- C, with C 1.39 
(see Eq. (|38p ). a change of the thermopower by tens of 
percents will cause a large change in the deduced value 
of X. These last fits should be regarded only as demon- 
strating our recommendations for a possible extension of 
the analysis of future experimental studies. 

At higher temperatures, the analysis will be influenced 
not only by data that are not in the critical region, but 
two further relevant physical processes may well come 
in. Obviously, inelastic scattering (by both phonons and 
other electrons) will be more important. Moreover, for 
T > ji — Erm some of the transport will occur via hopping 
of holes below the mobility edge. Their thermopower 
might cancel some of the contribution of electrons above 
/J, and thus reduce the thermopower below the values con- 



1.5 - 
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FIG. 6: A fit of the theoretical expression of Eq. (|26() to the 
experimental data. The data was taken from Ref. The 
thermopower is measured in units of fcs/e ~ 86^V/K. 



sidered here. This clearly needs further treatment. 

At any rate, the interactions appear to be strongly 
relevant and may give unexpected values for x. Mea- 
surements at lower temperatures and closer to the tran- 
sition are clearly needed. Checking simultaneously the 
behavior of both the conductivity and the thermopower 
is suggested as the method of choice for this problem. 

It must be mentioned that similar experiments on 
granular Al did not show the expected behavior. It 
should be kept in mind that the resistivities needed to 
approach the MIT for ^ lOOA grains are larger— than 
the ones for microscopic disorder. Another relevant issue, 
which we are going to examine in detail in future work, is 
that while Ep is in the eV range, all the energies (with- 
out electron-electron interactions) relevant for localiza- 
tion are smaller by several orders of magnitude than for 
microscopic disorder. Thus, the temperature range for 
the enhanced thermopower might well be in the sub K 
range, which was not addressed in the experiments. The 
Coulomb blockade may partially alleviate this, but only 
when it is operative (not in the metallic regime). 

The later experiment of Ref. [sl] , on Si : P obtained 
very modest enhancement of S, but were stated to have 
been dominated by effective magnetic impurities, which 
are known to be strongly relevant for the Anderson tran- 
sition (e.g. eliminating the weak localization contribu- 
tions). All these issues have to be clarified. 



V. CONCLUDING REMARKS 

In this paper we have shown how to include the ther- 
mal current in the Thouless scaling picture of conduction 
in disordered systems. Expressions were given for the 



1.4 




FIG. 7: The 4 data sets of Fig. [6] which are closest to the 
transition are shown, scaled to lie on a single universal curve, 
as function of T/Em- x was taken as 0.1 and each curve was 
given a single value of Em- 
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2x2 matrix of longitudinal thermoelectric coefScients, 
in terms of c7q{E), the T — conductivity of the system 
were its Fermi level fixed at the energy E. The Onsager 
relations were shown to hold within this formulation. For 
the usual critical behavior of (Jo{E), given by Eq. ([T]), 
these behaviors were analyzed for an arbitrary ratio of T 
to the distance to the mobility edge. They were shown 
to satisfy scaling relationships which were confirmed nu- 
merically along with their limiting behaviors. 

It was shown how the conductivity and thermopower 
data close to the Anderson QPT should be analyzed si- 
multaneously to yield a better estimate of the critical 
exponent x than the determination based on a(T) alone. 
This was done for the low-temperature limits of existing 
datSF^"^ on the transition in In20^-x, giving already a 
good ballpark estimate of x. The data going to higher 
(probably too high) temperatures do scale and collapse 
according to Eq. ([27l) . but the resulting values of x ap- 
pear to be too small. One may speculate that this is due 
to interaction effects, but we prefer to postpone this to 
after having done this analysis with lower temperature 
data closer to the transition. 

Similar experiments on granular Al do not produce a 
large and interesting thermopower as above. This is cer- 
tainly a matter for concern. The explanation might well 
be due to the smaller microscopic conductivity scales 
{^,R being the grain size), or to the different energy 
scales relevant for these systems <^ Alternatively, the in- 
elastic scattering, not treated in this paper, may be rel- 
evant as well. 

The sharp and asymmetric behavior of (Jo{E) near the 
transition is ideal for getting large thermopowers. The 
predicted values approach ~ 200^. While the experi- 
mental results-- on In20^-x are smaller by ^ 40%, this 
is still encouraging. Were it possible to increase these val- 
ues say by phonon drag effects, this might even become 
applicable. Clearly, a treatment of the effects of inelas- 
tic scattering on the thermopower is called for, especially 
including the hopping conductivity regime. 



where to get the last equality we used the equilibrium 
/ = (er ip 1)~^. The outgoing heat current TS is the 
time derivative of the population times the excitation 
energy. Thus, each particle leaving the system carries 
"on its back" an amount of heat E which is its energy 
(measured from /x). Summing Se over all energies shows 
that the outgoing heat current is given by the outgoing 
particle current where the contribution of each energy is 
multiplied hy E — ji. 

It should be noted that the equality of the amounts 
of E — ^ and —TS carried by the excitation implies 
that the relevant free energy does not change when the 
(quasi)particle moves to another system which is in equi- 
librium with the first one. This is true for equilib- 
rium fluctuations and also for linear response transport 
(y — > and AT — > 0), between the two systems. 

As remarked, the result that the heat carried by an 
electron is given by its energy measured from the chemi- 
cal potential, /i, is valid also for bosons. As a small appli- 
cation, one can easily calculate the net thermal current 
carried by a single-mode phonon/photon waveguide fed 
by thermal baths at T ± AT/2. The result is a thermal 
conductance of fc^T7r/(6/i) (per mode), with no reflec- 
tions. This agrees with the result of dl] . The sound/light 
velocity cancels between the excitation velocity and its 
(ID) DOS, exactly as in the electronic case. This is why 
this result and the one based on the Wiedemann-Franz 
law for electrons are of the same order of magnitude. 
That their numerical factors are equal is just by chance. 
With reflections due to disorder, once the waveguide's 
length is comparable to or larger than the localization 
length (mean free path for a single mode), its thermal 
conductance drops markedly. 



Appendix A: The heat carried by a transport 
quasiparticle 

To make this analysis useful also for heat transport by 
phonons, etc., we display the equations for both fermions 
and bosons. The entropy associated with a state of a 
given equilibrium system at energy E, having a popula- 
tion / is 

SE = -kB[flnf+{l±f)ln{l±f)], (Al) 

where the upper (lower) sign is for bosons (fermions). 
When the population f changes with time, the change of 
Se with time is 

Se = -ksfln-^ = -^/, (A2) 
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